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Abstract: 

We calculate the one-loop twist-3 gluon contribution to the flavor-singlet structure func- 
tion g 2 (x,Q 2 ) in polarized deep-inelastic scattering and find that it is dominated by the 
contribution of the three-gluon operator with the lowest anomalous dimension (for each 
moment N). The similar property was observed earlier for the nonsinglet distributions, 
although the reason is in our case different. The result is encouraging and suggests a sim- 
ple evolution pattern of g 2 {x, Q 2 ) in analogy with the conventional description of twist- 2 
parton distributions. 
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1. Introduction 



Twist-three parton distributions in the nucleon are attracting increasing interest as 
unique probes of quark-gluon correlations in hadrons. They have clear experimental sig- 
natures and give rise to various spin asymmetries in experiments with polarized beams 
and targets. Quantitative studies of such asymmetries are becoming possible with the 
increasing precision of experimental data at SLAC and RHIC, and can constitute an 
important part of the future spin physics program on high-luminosity accelerators like 
ELFE, eRHIC, etc. 

The structure function g 2 (x,Q 2 ) in polarized deep inelastic scattering presents the 
classical example of a twist-3 observable and received considerable attention in the past. 
The experimental studies at SLAC Jl|, 0, |3| have confirmed theoretical expectations about 
the shape of g 2 {x, Q 2 ) and provided first evidence on the most interesting twist-3 contri- 
bution. On the theoretical side, a lot of effort was invested to understand the physical 
interpretation of twist-3 distributions (see e.g. ||, ||, || for the review of various as- 
pects) and their scale dependence [0, [|, [LC], [11], [12|, |13[ . Conditions for the validity 
of the Burhardt-Cottingham (BC) sum rule J Q dx #2(2) = were discussed in detail 
151 , 0) Ol an d the second moment J^dxx 2 g 2 (x) was estimated using QCD sum rules 



18fl and on the lattice |19 



In spite of the significant progress that has been achieved, understanding of the scale 
dependence of g 2 (x,Q 2 ) still poses an outstanding theoretical problem. To explain the 
difficulty, we remind that to the tree-level accuracy the structure function g2(x,Q 2 ) or, 
equivalently, gr(x, Q 2 ) = gi(x, Q 2 ) + g 2 (x, Q 2 ) is given by the quark- ant iquark light-cone 
correlation function in a transversely polarized nucleon (see, e.g. |T5j) 



1 f°° d\ 

9t{x) = 2M J ^^^^^^^^(An)^^), (1.1) 

where p M and s M are nucleon momentum and spin vectors, respectively, and we assumed 
that the nucleon is moving in the z— direction,^ p M = (p + ,p^, 0±), p 2 = 2p + p^ = M 2 and 
p ■ s — 0. The light-like vector n is given by n M = (0 + , l/p+, 0±) so that n 2 = 0, (pn) = 1 
and the transverse direction is defined as orthogonal to the (p, n) plane. For comparison, 
the leading twist-2 spin structure function is written as 

/°° dX 
— e iXx (p, s z \m il^n)\p, s z ) . (1.2) 

In the parton model, g\ (x) measures the quark helicity distribution in a longitudinally 
polarized nucleon. Such an interpretation can be made because the quark helicity opera- 
tor E p commutes with the free-quark Dirac Hamiltonian H = a z p z . On the contrary, the 
quark spin operator projected along the transverse direction does not commute with 
the Hamiltonian and thus there exists no energy eigenstate \p z ) such that Ej_|p 2 ) = s±\p z ). 
The transverse spin of the nucleon cannot, therefore, be thought of as being composed of 
transverse spins of individual quark (gluon) constituents. The transverse spin average of 

1 Throughout the paper we shall use the following definition of the light-cone components p± = 
(po ±p3)/y/2 and p± = (pi, Pi)- In addition, we shall not display the gauge factors connecting the quark 
fields and ensuring the gauge invariance of nonlocal light-cone operators. 
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quarks in ( |1 . 1| ) that defines gr{x) is sensitive to the dynamics of quark-gluon interactions 
and does not have any probabilistic interpretation in terms of one-particle quark parton 
densities. 

One possible way to see the relation of the transverse spin to gluonic degrees of freedom 
is to decompose the quark field operator in "good" (+) and "bad" (-) components q(x) = 
q+(x) + q-(x) where P + = ^7-7+ and P_ = ^7+7- are the corresponding projection 
operators ||20|| . It is easy to check that the correlation function in ( |1.2| ) involves only 



good quark components, while in ( |1 . 1|) necessarily one good and one bad components 
are involved. In the approach of [^(| only good field components correspond to genuine 
partonic degrees of freedom, while bad components are not dynamically independent and 
can be eliminated through the equations of motion in favor of good components and 
insertions of quark masses or gluon fields. An important point is that this relation is 
nonlocal and involves quark and gluon fields with different positions on the light-cone 



ff-(An) = J^e-^-^^ ± (X'n) q+ (X'n), (1.3) 

where = <9 M — igA^ is the covariant derivative. 

Eq. ( |1.3| ) states that degrees of freedom associated with the bad component of the 
quark field in ( |1.1| ) are, in fact, those of one quark and one gluon. The structure func- 
tion gr{x, Q 2 ) is, therefore, naturally related to the quark-antiquark-gluon correlation 
function in the nucleon. More precisely, gr{x) presents by itself only one special projec- 
tion of this more general three-particle distribution -D(£i,£2,£s) that depends, generally, 
on the momentum fractions £$ carried by three partons. It is this special projection, 
grix), that can be measured in deep inelastic scattering with a transversely polarized 
target. On the other hand, the scale dependence of the quark-antiquark-gluon distri- 



bution function involves the "full" function Z?(£i, £2^3) in a nontrivial way || |2l], ^] 
and the knowledge of one particular projection gx(x, Qq) at a given value of Ql does not 
allow to predict gr{x, Q 2 ) at different momentum transfers: a DGLAP-type evolution 
equation for gr(x, Q 2 ) in QCD does not exist or, at least, is not warranted. The reason 
is simply that inclusive measurements in general do not provide complete information on 
the relevant three-particle parton correlation function -D(£i, £2, £3)- 

From the phenomenological point of view such situation is not satisfactory since it 
would mean that one cannot relate results of the measurements of g2{x) at different values 
of Q 2 to one another without model assumptions. The theoretical challenge is, therefore, 
to find out whether the complicated pattern of quark-gluon correlations can be reduced 
to a few effective degrees of freedom. One has to look for meaningful approximations to 
the scale dependence that introduce a minimum amount of nonperturbative parameters. 

In particular, it was found by Ali, Braun and Hiller (ABH) |H| that the scale- 
dependence of the flavor- nonsinglet contribution to gr(x, Q 2 ) simplifies dramatically in 
the limit of large number of colors N c — > 00. To explain this result, it is convenient to 
use the language of the Operator Product Expansion (OPE), see Sect. 2 for more details. 
The statement of the OPE is that odd moments n = 3, 5, . . . of the structure function 
#2 (x, Q 2 ) can be expanded in contributions of multiplicatively renormalized local twist-3 
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quark- ant iquark-gluon operators^ 



„1 n-3 , 

/ dxx n - 1 g 2 (x,Q 2 ) = J2cLA 



^y-'aoM)), (i.4) 



where C^_ 3 are the coefficients and ((0^_ 3 (/x))) reduced matrix elements normalized at 
the scale \i\ 7^_ 3 are the corresponding anomalous dimensions that we assume are ordered 
with k: 7°_ 3 < 7*_ 3 < . . . < 7^l| for each n, and 6 = lliV c /3 — 2rij/3. The exact analytic 
expression for the lowest anomalous dimension in the spectrum, 7°_ 3 , has been found in 
Til an d it was also noticed that the coefficient functions of all other operators (with 



higher anomalous dimensions for each n) are suppressed by powers of 1/iV 2 . Thus, to 
the stated 0(1/ N^) accuracy, each moment of #2(2) involves a single nonperturbative 
parameter while the complicated degrees of freedom related to quark-antiquark-gluon 
correlations essentially decouple. The result can be reformulated as a DGLAP-type 
evolution equation 



4Cpl + S(1- Z 



P NS (z) 



1 / 7T 



2C F , (1.5) 



1 - z_ 

where Cf = (JV| — 1)/(2N C ) and we have included the 1/N% corrections calculated in 



The present paper is devoted to the extension of this analysis to the flavor-singlet 
sector in which case twist-3 three-gluon operators have to be included. We calculate the 
leading one-loop 0(a s ) gluon contribution to the coefficient function and examine its 
properties. We find that the one-loop coefficient function is such that it mainly picks up 
the contribution of the twist-3 three-gluon operator with the lowest anomalous dimension 
for each moment N. The dominance of the lowest three-gluon "state" is observed both for 
the logarithmic contribution ~ lnQ 2 //i 2 that reflects mixing with the quark-antiquark- 
gluon operators, and the constant term that gives rise to a "genuine" gluon contribution. 

The result is very encouraging and allows to hope for the similar pattern of a simplified 
evolution that mainly involves a single quark-antiquark gluon and a single three-gluon 
parton distribution corresponding to the "trajectories" with the lowest anomalous di- 
mension as important degrees of freedom, although the reason for such a simplification 
is different. The properties of these trajectories have been recently studied in |2~T| . 

The outline of the paper is as follows. Sect. 2 is introductory and reviews existing 
results on the OPE of the antisymmetric part of the T-product of two electromagnetic 
currents to twist-3 accuracy. The calculation of the coefficient function of twist-3 three- 
gluon operators, is presented in Sect. 3 and its structure is elaborated upon in Sect. 4; 
in Sect. 5 we summarize. Technical details on the twist separation in gluon operators are 
presented in the Appendix. 



2 Here we neglect the twist- 2 contribution to (^(x). 



2. The Operator Product Expansion 

As well known, the hadronic tensor which appears in the description of deep inelastic 
scattering of polarized leptons on polarized nucleons, involves two structure functions 

1 



W {A) 



p-q 



s p p p 

p-q 



92(xb,Q 



(2.1) 



where the nucleon spin vector is defined as = u(p, s)7 M 7 5 -u(p, s), u(p, s) is the nucleon 
spinor (u(p, s)u(p, s) = 2M, s 2 = — 4M 2 ), and is related to the imaginary part of the 
Fourier-transform of the T-product of two electromagnetic currents, antisymmetrized 
over the Lorentz indices: 



iTfi = \J d'xe^(p,s\T{ 3ll (x/2)j u (-x/2) - j u (x/2) 3tl (-x/2)}\p,s} . (2.2) 

We are going to examine the light-cone expansion of fl2.2|) at x 2 — > and write down the 
answer in terms of nonlocal light-cone operators of increasing twist, schematically 

^ T {j^( x )ju(-x) -j u (x)j lx (-x)} x = ° 

'-J^^-JIC * 0,r- 2 + [C * 0,] tw - 3 + (higher twists)} , (2.3) 

where C * Op stands for the product (convolution) of the coefficient functions and oper- 
ators of the corresponding twist. This expression is explicitly [/(l)-gauge invariant, i.e. 
d/dx,T{j,(x) 3u (-x)} = 0. 



2.1. Leading-order results 

To the leading order in the strong coupling, the OPE of the antisymmetric part of the 
T-product in ( [2.3| ) can be written in a compact form as [2^ 



\c*o,r- 2 = x ^ y, 



q=u,d,s, 



du q{ux) fysq^—ux) + (x — > —x) 



i i>u 

du j dv [(it + v) St(u, v, —it) + (u — v)Sp(i 



q=u,d,s,... 

+ (x — > —x) , 

where notation was introduced for light-cone nonlocal quark-gluon operators 

i±r„ i. ^ _ X ^„„\\A„n tu„\^„n ti^s„^j„t^.\ (2.5) 



0] 

(2.4) 



(a, b, c) = -q(ax)[igG flu (bx) ± gG ^ v (bx)~f b ]x v tfq(cx) . 



The dual gluon strength tensor is defined as G^ u = \e ^a^G 0113 and we use the conventions 
75 = ^7o7i7273 and e 0123 = 1, see p4| . To save space, here and below we do not show the 
gauge factors connecting the quark (gluon) fields: 



[ax, bx] = P exp 



ig / dux^A^(ux) 



(2.6) 



a 



b 



Figure 1: Leading-order twist-3 contributions to the OPE of the T-product of electromagnetic 
currents 

The full twist-3 contribution ( |2.4j) to ( |2.3| ) is assembled from the Feynman diagrams 
shown in Fig. [I]. The individual contributions of the two diagrams in Fig. [I] correspond 
to the two possibilities to apply the derivative in (|2.3| ): d / dx a {S^ / x 2 ) = —2x a S^/x 4: + 
l/x 2 dSg /dx a . The term ~ I fx 4 presents the contribution of the the diagram in Fig. |l]a, 
rewritten in terms of quark-gluon operators using equations of motion. The contribution 
~ 1/x 2 corresponds to the diagram with gluon emission from the hard propagator in 
Fig. [TJb and is necessitated by gauge invariance of (|2.3|). 

Going over to matrix elements, we introduce the usual quark helicity distributions 
Aq(x B ) = q\x B ) - q l ( x s) 

(p,s\q(x) fr 5 q(-x)\p,s) = (sx) J d£ e 2i ^ px Aq(£, p, 2 ) , (2.7) 

where \x refers to the normalization scale of the operator in the l.h.s. and positive (nega- 
tive) £ correspond to the contribution of quarks (antiquarks) Aq(x B ) = q^(x B ) — q^(x B ), 
Aq(—x B ) = Aq(x B ) = q\x B ) — ^(ib), respectively. 

Similarly, we define the twist-3 quark-antiquark-gluon parton correlation functions as 
(cf. ' 

(p,s\S±(u,v,-u)\p,s) = 2i{px) [s^x) -p^sx)\J VZe^^-^Dfi^^^) 

(2.8) 

where ^ 

J Vi = J d£id£ 2 d£ 3 5(6 + 6 + 6) • (2.9) 

The correlation functions _D^(£j) have the following symmetry property: 

^(6, £2, 6) = P g T )*(-6,-e 2 ,-ei) (2.10) 

They are in general complex functions, but the imaginary parts do not contribute to the 
structure functions and can be omitted . 

3 We always imply that x 2 can be put to zero in the operator matrix elements. Under this condition, 
the nonlocal operators Si still contain a superposition of twist-3 and twist-4 terms [ p6| , with the twist-4 
terms being explicitly proportional to Xp. The easiest way to separate the genuine twist-3 contribution 
is to take the transverse projection — ► S±, see also || Q for explicit expressions. 



For further use, it is convenient to introduce a separate notation for the nucleon 
matrix element of the specific combination of quark-antiquark-gluon operators entering 
the OPE in (|2l 



dv (p, s\ [(u + v) S+(u } v, -u) + (u- v)£ (u, v, -u)] \p, s) 



(sxj 



(px) 



J d( . e mu px A g T (^ 2 ). (2.11) 



The function AqT{xB,Q 2 ) will play an important role in what follows. It describes 
the momentum fraction distribution inside the nucleon of the transverse spin and has 
the same support property as the parton distribution in (|2.7|) . However, in contrast 



with (|2.7| ), it does not have any probabilistic interpretation but can rather be expressed 
through the more general three parton correlation functions _D^(£i, £2, £3), integrating out 
the dependence on the gluon momentum fraction ||. Making a Fourier transformation 
of (|2.4| ), taking imaginary part and comparing with the definition of structure functions 



in ( |2.1|) one obtains 

9l (x B , Q 2 ) = \Y, 4 t A ^> f = Q 2 ) + M-*b, /i 2 = Q 2 )} , (2.12) 



(2.13) 



(2.14) 



where xb = Q 2 /(2pg) is the Bjorken variable, and |23|, [H| 

92(x B ,Q 2 ) = g? w (x B ,Q 2 ) + f -l^T(y,Q 2 ) + Aq T (-y,Q 2 )} 

A q Jx B y 

where 

g? w (x B , Q 2 ) = - 9l (x B , Q 2 ) + f 1 ^ 9l (y, Q 2 ) 

Jxb y 

is the the familiar Wandzura-Wilczek contribution |27J and will not be considered from 
now on. 

Going over to the moments in Q2.13|) we obtain 

/ dx B x B - 1 g 2 (x B ,Q 2 ) = ^-Y^e 2 q [ dx B x^,~ 1 [Ag T (xB, Q 2 ) + Aq T (—x B , Q 2 )} . (2.15) 
Jo Zn Jo 

For odd n > 3 the relevant integrals of the transverse spin quark distributions Aqx 
are given by the reduced matrix elements of local twist-3 operators that arise via the 
Tailor-expansion of the nonlocal operators ( |2.5| ) at short distances x M — > 0: 

[S%\n = \l(D- xf[igG, u ± gG^js] j;x v (5- x) N ~ k q . (2.16) 
According to ( p.8|) , the reduced matrix elements ((. . .)) of these operators 



MTOiP,*) = 2(ipx)» +l Mpx)-p,(sx)} «[$*]*,» (2.17) 
are equal to moments of the quark-antiquark-gluon three-particle distribution amplitude 



J^^t k Dt(Ci). (2.18) 



The symmetry relation ( |2.10| ) implies 



(2.19) 



Expanding Eq. ( |2 . 1 1 ) at short distances we obtain 

-1 N 

/ dU N+2 Aq T (0 = 2Y,(-l) N - k [(k + miS-1%)) + (N-k + l)(([S+] k N ) 

J —1 i — n 



k=0 

N 

k=0 



\N-k 



4 I VtReD-{t; <} 



1 



<9£i 6 + 6 



JV+2 



\AM-2 



(2.20) 



The last equality can also be obtained directly from the definition in ( |2.8| ), (|2.11 ). The 
following comments are in order. 

We note the function Agy(£) takes real values. The expression in the last line of 
( p.20| ) can be used for an analitic continuation to N — > —2 and remains finite provided 
that the corresponding integral of the D— function converges. This convergence, thus, 
presents a necessary condition for the validity of the BC sum rule at N = —2. 

At N = — 1 an absence of a local twist-3 operator ( |2.16 ) with dimension four implies 
the constraint 

J dZZAqr(0 = J dx B x B [Aq T (x B )-Aq T (-x B )] = 0. (2.21) 

This relation should be compared with the first moment of g2{x) given by (p. 15 ) that 
involves the combination of the same distributions but with a different C-parity 

/ dx B x B g 2 (x B ,Q 2 ) = - y^el f dx B x B [Aq T (x B ) + Aq T (-x B )] ■ (2.22) 
Jo 4 _ 

Vanishing of this integral (known as Efremov-Leader-Teryaev sum rule |28| , ^|) is, there- 
fore, not warranted by the OPE, although its numerical value can be small since the r.h.s. 
does not receive contribution from the valence quarks. 



2.2. The scale dependence 

The dependence of the structure function g2(x,Q 2 ) on Q 2 is driven by the scale de- 
pendence of the distribution functions Aq T (x, n 2 ). Going over to moments ( |2.2(J| ), it 
corresponds to the renormalization-group scale dependence of the local operators [S^]^. 
Similar to the familiar case of the helicity distributions Aq(x B ), one has to distinguish 
between the components with different flavor symmetry as they have a different scaling 
behavior. For instance, the flavor decomposition of the u-quark distribution looks like 

Au T {x B ) = I (Au T - Ad T ){x B ) + I {Au T + Ad T - 2As T ){x B ) 
2 o 

+ - (Au T + Ad T + As T ){x B ) . (2.23) 
o 



Renormalization of flavor- nonsinglet contributions, Ag NS , given by either Au — Ad, or 
An + Ad — 2As is simpler since they do not mix with gluons. Still, the number of 
contributing operators rises linearly with iV 0. By an explicit calculation one obtains, 
for the two lowest moments @, ||, [H| |ll| 



j/Ve 2 Agf(e,g 2 ) = ^ 0/6 (^o (^ 2 ))) 



dte^ s (t,Q 2 ) = L^ b 



(2.24) 



0.807<(S 2 V)» - 2.320«S 2 V)» + 2.894«S 2 V)» 
+ L^ .028((S 2 V)>> -0.014«S*(/i a )» - 0.026«5*(/* 2 )» 
+ L^Ho.l65<<S 2 V)>> + 0.334«S 2 V)» + 0.132«£ 2 V)» 



where we used {{S^(fi 2 ))) = Re (([5 ]%(^ 2 ))) as a shorthand and L = a s (Q 2 )/a s (n 2 ). 
Anomalous dimensions are equal to |10[ 

(2.25) 



7o ° = 8.5, 72 = 10.89, T £ = 13.71 , 7 2 = 16.15. 



Eq. (|2~24p illustrates the main difficulty: since the number of contributing operators 
prolifiterates with N, so does the number of independent nonperturbative parameters 
((S^(fi 2 ))). On the other hand, Eq. (|2.24|) reveals a remarkable pattern: coefficient func- 



tions in front of the operators with higher anomalous dimensions are much smaller than 
those with the lowest anomalous dimension. This structure is not accidental, but related 
to a dramatic simplification of the renormalization-group evolution of flavor-nonsinglet 
operators in the large— N c limit. As was found in [101, the small coefficients in (|2.24|) are 
in fact suppressed by powers of 1/iV 2 and one obtains in the limit N c — > oo 



dtt 2 Agf (£, Q 2 ) N °=°° L^ S / b \((S^ 2 ))) - 2«S 2 V)» + 3«S 2 V)» 



(2.26) 



so that the scale evolution of the second moment of Aq^ involves the same moment. 
The similar phenomenon takes place for arbitrary N: 



d^ N Aq^,Q 2 ) N °=°° L^J^^Aq^,^) . (2.27) 



Here 7^ s is the lowest anomalous dimension in the spectrum of flavor-nonsinglet twist-3 
operators. It is known analytically in the large— N c limit [[K|] and the 1/N 2 corrections 



have been recently calculated using the large— iV expansion in |LJ 

1 



N c [2^{N + 3) + 2 lE + 



N + 3 



2 'ln(iV + 3)+ 7E + l-^ + (l/iV 2 )) + O (Jp) . (2.28) 



EOM 




a 



Figure 2: Leading-order twist-3 gluon contribution to the OPE of the T-product of electro- 
magnetic currents. 



Here ip(x) = d\nT(x)/dx stands for the Euler ^-function. Eq. (|2.27|) together with (2.28 ) 
are equivalent to the DGLAP evolution equation in (|1.5| ). 

The scale dependence of the flavor- singlet distribution Aqf. = Aut + Ag?t + Ast 
differs from the above in three aspects. First, the mixing matrices of the relevant quark- 
ant iquark-gluon operators receive extra terms related to the possibility of quark-antiquark 
annihilation. Second, they mix in addition with an entirely new and equally big set of 
three-gluon operators. Third, the three-gluon operators themselves contribute to the 
OPE of the T-product of the electromagnetic currents starting order a s . Concerns have 
been raised (see e.g. [p9fl ) that in particular the last contribution does not have a simple 
structure and will spoil any ABH-type approximation in the singlet case. We begin, 
therefore, with the corresponding calculation. 



3. The gluon contribution to the structure functions 



The leading-order gluon contribution to the T-product of two electromagnetic currents in 
( j2.2p is described by the box diagram shown in Fig. |2|. Its calculation can be easily done 
using the background field approach of Ref . || . Namely, considering outgoing gluons as 
classical background Yang-Mills fields we calculate the box diagram replacing the free 
quark propagators by propagators in an external fields. Then, the antisymmetric part of 
the T-product in Q2.2|) is given by 



H$\x, -x) = -T{ Jfl (x)j u (-x)} -(p~v) 



9 Yl e S Tr J»S{x,-xhuS{-x,x) -{n^v) 



(3.1; 



q=u,d,s, 



where Tr denotes the trace over both color and spinor indices, and S(x, —x) is the quark 
propagator in a background gluon field 



S(x, —x) = (x\ 



— x) . 



(3.2) 



For x 2 — > the propagator S(x, —x) exhibits light-cone singularities that one handles 
using the dimensional regularization with d = 4 — 2e and e < 0. Then, expanding 
S(x, —x) in powers of the deviation from the light-cone and retaining contributions of 
gluon operators up to twist-3 we find 



(4ii) d / 2 S{x, -x) 



T{d/2) T(d/2 - 1) r(d/2 - 2) 



First two terms in the light-cone expansion of the propagator are known 



(3.3) 



Si(x, -x) = - J du{ [x,ux] gG xa (ux)-f a -f 5 - iugG xa (ux)^ a 



[ux, —x\ 



+ J dv (1 —u)(l + v)[x,ux]gG xa (ux)[ux,vx)gG xa (vx)[vx, -x]j (3.4) 

and the further ones can be calculated using the technique described in Appendix A 
of ||. Here and below we use a shorthand notation G xa = G a ^ a x^t a with t a being the 
generators of fundamental (quark) representation of the SU(N C ). 

Omitting the disconnected (gluon field independent) contribution ~ l/(—x 2 ) d , one 
obtains 



n<£(*,-*) = 



T 2 (d/2 



q=u,d,. 



' 9 (47r) d (-x 2 ) d - 2 



■^■/j,i>\.Xj X) ~\~ j ^B^v^x, x) 



-0 



-X 



2 \d-3 I ' 

(3.5) 



where 



(3.6) 



l t ,S 2 {x, -x)y v S (-x,x) 



Tr 



7 At S (a;, -x)j u S 2 (-x,x) J - (ji v) 



Several comments are in order. First, the expression in (|3.5|) defines the most singular, 
~ 1/x 4 as d — > 4, contribution to the light-cone expansion of the T-product (|2 . 3| ) . This 
contribution alone suffices to determine the coefficient functions of the twist-2 and twist- 
3 gluon operators since less singular contributions of the same twists can be uniquely 
restored from matching (|3.5| ) into the general ?7(l)-gauge invariant expression (|2.3| ), cf. 
the discussion in Sect. 2.1. Second, notice that the first term in ( |3.5|) is analytic and the 
second is singular in the limit d — > 4. These two terms correspond to the two distinct 
integration regions in the quark momentum in the loop (see Fig. 2) k 2 ~ l/(— x 2 ) ~ Q 2 
and fc 2 < Q 2 , respectively. The first term, coming from large momenta, determines the 
one-loop 0(a s ) coefficient function of gluon operators at a hard scale of order Q, while 



the second term will be interpreted as a tree-level quark coefficient function times the 
one-loop evolution (mixing) into gluons. Finally, one can convince oneself that the traces 
in (|3.5|) can be calculated in dimension d — 4. This is obvious for the first term, and can 
be shown for the second, the reason being that calculation of diagrams of the type shown 
in Fig. 2 does not involve contraction of Lorentz indices of 7-matrices. 

Calculating ( |3.5| ) we shall assume the translation invariance of U.^ A '(x, — x) along 
the light-cone[J In addition, we shall impose the equations of motion for gluon fields, 
[D^,G^ u (x)} = 0, which amounts to putting external gluons in the box diagram (see 
Fig. 2) on their mass-shell. Then, using the explicit expression for Si(x, —x) given in 
ni) one obtains 



A^vix, —x) = 8g x^l duuuTr <G x ^(ux)G u ^(—ux) — G x ^(—ux)G^(ux) 



+ Aig Xnl duu 



dv 



(l+v-^J O v (u, -u,v) + (l-v-^j O v (v,u, -u) 

(3.7) 

where we used u = 1 — u etc. and introduced a notation 

O v (u,v,t) = gTr | G xa (ux),G xu {vx) G xa (tx) 



+ -uO u {u,v,-u) -(fi^u) 



'jf abc G a xa (ux)Gl(vx)GUtx) 



(3.8) 

In the last relation the gauge factors between the gluon fields in the adjoint representation 
are implied. 

The calculation of the second, singular contribution in (|3.5| ) is more tedious. After 
considerable algebra we obtain, however, an equally simple expression 



Aig x^l duu I dv I 3+v — u O u (u, -u,v) + 3-v — u J O v (v , w, — u) 



+ -(2 + u) O v {u, v, -u) -(n^v). 



(3.9) 



Substituting Eqs. ( |3.8| ) and ( |3.9| ) into ( p.5| ) we finally obtain the leading order expression 
for the antisymmetric part of the T-product of electromagnetic currents that takes into 
account both twist-2 and twist-3 contributions. However, in order to match ( |3.6|) into 
the general OPE form ( [2.3|) we have to separate the different twists. 



3.1. Twist separation 



The two-gluon operators in ( |3.71 ) and ( p.9[ ) contain, generically, contributions of both 
twist-2 and twist-3. The separation of twists corresponds in this case to the separation 
of the terms of different symmetry and can be done using the trick described in || . The 



4 That is, neglect contributions of the operators containing total derivatives that have vanishing for- 
ward matrix elements. 



twist-2 part of the two-gluon operator in ( |3.7| ) can be written as (see Appendix A for 
details) 



Tr {G xi {x)G^(-x)} 



tw-2 



d 

duu— — Tr {G x ^{ux)G x ^{— ux)} . (3.10) 
o ox v 



Similarly, neglecting the irrelevant operators proportional to the equations of motion and 
total derivatives, one obtains the twist-3 contribution as 



Tr {G x ^(x)G^(-x)} 



tw— 3 



*1 pu 

il du I dv Tr <| (1 + u 2 ) [G x ^{ux),G x ^{vx)\G xvK 



-ux) 



-(1 - uv) [G x z(ux),G xl/ (vx)]G x t(-ux) 
+ (1 - u 2 ) [G xv {ux),G x z(vx)]G x6 {-ux) }. (3.11) 



Applying ( |3.10|) we find that the twist-2 two-gluon contribution to (|3.7| ) in fact cancels 



out, and the remaining twist-3 part can be rewritten using ( |3.11| ) as 
[Tr {G x ^ux)G^{-ux) - G x ^-ux)G u!i (ux)}}^ = 

= -2i J ds J dt{o u (s,-s,t) +O u (s,t,-s) +O u (t,s,-s)\. 



(3.12) 



Substitution of this relation into (|3.7| ) yields 



[A 



ltw-2 







/*1 r-u 

[A^ u ] tw ~ 3 = —Aig 2 x fl duu 2 dv \ (u — v) O u (u, — u, v) + 2u O v [u, v , — u) 

JO J-u 



+ (U + V) O u (v, U, —U) \ - (yU <->• v) . 



(3.13) 



The separation of twists in (3J3) is equally simple. To this end we note that (in dimension 
d = 4) 

G x Jux)G xu (-ux) - G xu (ux)G x J-ux) = e^ upa x a G xa (ux) G pa (-ux) (3.14) 



and the relations ( |3.10|) - (|3.11| ) remain true to the claimed accuracy if one of the gluon 
strength-tensors is substituted by its dual counterpart. We obtain 

[B^] tw ~ 2 = -16g 2 e f , up(7 x (T -^— \ du (ulnu + uu)Tr {G x ^(ux)G x ^(-ux)} , 



[B, 



ltw— 3 



I6ig 2 x^j' duj dv j(uw + ju 2 + u\nu) v O u (v, u, — u) + u O v (u, v , — u 



— v O u (u, —u, v) 



-—u 2 (u + 2) O u (u,-u,v) + O u (u,v,-u) + O u (v,u,-u) j 

(3.15) 



3.2. Twist-2: Results 



Substituting ( |3.15 ), ( |3.13| ) into ( |3.5|) , subtracting the (collinear) singularities in the MS 
scheme^ and combining with the leading-order result in ( |2.4j ) , we obtain the twist-2 
contribution 



[C * O p ] 



tw— 2 



X/3 s—^ 
X 1 ^ 



q=u,d,s, 



du 



q{ux) ^7 5 g(— mx) + (x — > —x) 



MS 



(3.16) 



4a, 



TT 



H — ( ln(— x 2 fiL_) + 2^e ) (win it + uu)Ti {G x ^(ux)G x ^(—ux)} 



where the subscript [. . .J^a indicates the normalization scale of the operator. Note sim- 
plicity of the answer: the entire gluon contribution can be eliminated by choosing the 
proper scale of the quark operator /i 2 -^ = l/(— x 2 e 27E ). This property is lost in the mo- 
mentum space since after the Fourier transformation contributions of different light-cone 
separations get mixed. 

Going over to the matrix elements, we introduce the usual gluon helicity distribu- 



tion EOL 31 



(p, s|Tr <^ G xa {x)G xc 



-x 



4 



sx)(px) / ci£e 2 ^A^,/i 2 ). 



(3.17) 



Note that Ag(0 = Ag(-£). 

Moments of the structure functions ( |2.1|) are obtained by the expansion of the T- 
product Q2.2|) in momentum space in powers of uj = —2(pq)/q 2 , Q 2 = —q 2 in the unphys- 
ical region to — > 0, and matching to the corresponding expansion in terms of structure 
functions: 



T$ = _4e 



- jlV 



uvaB r. 



n=l,3,.. 



CO 



n-l 



s (3 g 1 (n,Q 2 ) + 



s p pr 

p-q _ 



92(,n,Q 2 



where the moments are defined as (for any function /) 

f(n,Q 2 )= I dxBxl' 1 f(x B ,Q 2 ) 
Jo 



(3.18) 



(3.19) 



Taking the Fourier transform of (|2.3|) by using ( |3.16D and matching the obtained 
expression into (|3.18|) we obtain after some algebra 



9i(n, Q 2 ) = - Yl e l \ Aq ( n > ^Ms) + A ^( n ' ^Ms) 



q=u,d,s 



ft. 



n — 1 



[92(n,Q 2 )} 



2Mtw-2 



2ixn{n + 1) 
n — 1 



Ag(n, /i 2 



In 



01 



- ip(n) - 1 - 7£ 



J2_ 



gi(n,Q z 



(3.20) 



5 Note that the (non-renormalized) tree-level coefficient function in coordinate space contains a 
notrivial d— dependence corresponding to the quark propagator T(d/2) jt/{—A-n:x 2 ) d t 2 . Making the 
MS— subtraction in coordinate space one has to keep this factor in dimension d in the counter-term. 
Alternatively, one can make the Fourier transform first, and then subtract the divergencies in the usual 
way. One can check that the renormalized coefficient functions obtained in both ways are indeed related 
to each other by a 4-dimensional Fourier transformation. 



in accord with the Wandzura-Wilczek relation ( |2.14| ). Comparing the first expression in 
( p.20| ) with the general expression 



q=u,d,t 



9i(n,Q 2 ) = ^- e 2 q lAq(n,Li 2 ) + Aq(n,fi 2 ) + ^ qg (n) In — + C g {n) Ag(n; fi 



we find the anomalous dimension and the gluon coefficient function as 

a, n — 1 



7wW 
C g (n) 



2tt n(n + 1) ' 

a s n — 1 
27rn(n + 1) 



ij}(n) + 1 + 7s 



Going over from the moments to the momentum fraction representation, 
lqg( n ) = Jo dx x n ~ 1 P gg (x) and C g {n) = dx x n ~ 1 C g (x), we get 



C g (x) 



2^ 



(2x - 1) In 



x 



1 — x 



+ Ax - 3 



(3.21) 



(3.22) 



(3.23) 



Expressions in Q3.22 ) and ( 3.23|) are in agreement with the well-known results, see e.g. 



3.3. Twist-3: Results 



The twist-3 contribution to the T-product ( |3.5| ) comes from ( |3.13| ) and ( |3.15| ). To cast it 
into the U(l)— gauge invariant form ( |2.3| ) we notice that 

(3.24) 



Xfj,O u (u, v, t) - x u Of,(u, v, t) = E^ OL0 x a O g {u, v, t) + 0(x 2 ) , 
where the operators O v were defined in (|3.8| ) and 

Op(u,v,t) = ^r bc G a xa {ux)G%{vx)Gl a {tx) . 



(3.25) 



Then, combining together Eqs. fl3.13|) and (|3.15|) , subtracting the collinear singularities 
in the MS scheme (see previous footnote) and comparing with the general structure of 
the OPE in (O) we obtain: 



[C * O p ] 
i 



tw— 3 



V e 2 

2x 2 c i 

q=u,d,s,... 



duj dv | (u + v) Sp(u, v, —u) + (u — v)Sp(—u, v, u) 
+u 2 — \ (u + v)Op(v, u, —u) + 2uOp{u ) v, —u) + (u — v)Op(u, —u, v) 

7T L 



(3.26) 



I — — (ln(-x 2 fi 2 ^) + + l) (mm + \u 2 + ulnu) vOp(v,u,—u) +uOp(u,v,—u) 
1 



-v Op(u, -u, v) - —u (u + 2) 



Op(U, -U, V) + Og{U, V, -u) + Og(v, U, -11 



'a 



3[ 



1 1 



where the subscript indicates the normalization point of nonlocal quark and gluon 
operators. Here we introduced the C-even quark-gluon operator 

S^u, v, -u) = S+(u, v, -u) + S~(-u, v, u) (3.27) 

so that 

S^-u, v, u) = S~(u, v, -u) + S+(-u, v, u) . (3.28) 

The following comments are in order. 

The twist-3 gluon contribution in (|3.26|) has two parts. The lengthy expression in the 
last two lines can in fact be eliminated by the choice of scale in the quark operator in 
the first line: //|^ = l/(— x 2 e 2 ~ tE+1 ). As a nontrivial check of our calculation, we have 
verified that our answer ( |3.26|) is in agreement with the renormalization group equation 
for the twist-3 operator S M (w, v, —u) 0, [T2[ 



2 ru r-s 

[Sp(u 1 v,-u)] l * = [Sp(u,v,-u)] IJ 2 -— ^-ln-f/ ds dt(2uY' i 

^ ^\J-u J-u 

x{[2u(a -t) + 4(w - s)(t + s)}Op(s,v,t) - 2\u\(s - t)[Op(s, t, v) - Op(v, s, *)]}.(3.29) 

In addition, the expression in the second line in ( |3.26|) defines a 'genuine' twist-3 gluon 
coefficient function that cannot be eliminated by the scale choice in the quark operator. 
This expression is surprisingly simple and can be cast in the form similar to that of the 
tree-level contribution of the quark operators, Eq. (|2.4f) : 



(3.30) 



q=u,d. 

where 



Qp(u, v, -u) = -— u 2 Op(v,u,-u) + Op(u,v,-u) . (3.31) 



7T 



The nucleon matrix element of (|3.26|) defines the twist-3 gluon correlation function 
similar to (|2.8l )fl 

(p, s|0>, v, -u) \p, s) = -2{px) 2 [s^px) - p^sx)} j Vi e^^-^D^, &) 

(3.32) 

with the integration measure given by (|2.9| ). Since 0^(u,v, —u) = —0^(—u,v,u) ac- 
cording to the definition ( p. 251 ), the correlation function D g {l~i) is antisymmetric to the 
interchange of the first and the third argument: 

60 = -D g fa, 6, 6)- (3-33) 
Expanding the nonlocal operator ( |3.25| ) over local twist-3 gluon operators 



[G,) k N = ^9f abc G a xa (D- x) k G^ (D- x) N - l - k G c xa (3.34) 



3 cf. the footnote to (12.8 



i r 



and denning the reduced matrix elements ((...)) of these operators 

MP/JkM = 2(jpx) N+1 [a M (px) -P,(sx)] (([G] k N )) (3.35) 
we obtain the moments of the gluon three-particle distribution amplitude fl3.32|) as 



«m)) = J VZZitr Dgfo) (3-36) 

with < k < [N/2] - 1. The symmetry flOg ) implies that (([G] k N )) = - {{[G]%- 1 ~ k )) and, 
therefore, the number of independent gluon matrix elements is equal to [iV/2]. 

Finally, substituting (|3.26| ) into ( |2.3|) and taking the Fourier transform ( |2.2j ), we match 
the result into the expansion (|3.18|) to obtain the moments of the structure function for 
n = 1, 3, . . . 



q=u,d,s,... ^ 



) 



(3.37) 



+ 



o:. 



MS' 



47r n + 1 



<^(6)+^(6) ^ 



Q 2 



MS 



where, as usual, applicability to the lowest moments relies on the extra assumption 
about the high-energy asymptotics of the cross section, alias the assumption that the 
corresponding dispersion relation does not involve extra subtraction constants. Here 
-D 9 (6) is the distribution function corresponding to the C— even quark-gluon operator 
0) 

D q(&, f4ss) = Re D i (& ^Ms) ' (3-38) 
the quark coefficient function is defined as 

d d 



$£(6,6 



<96 



6 + 6 



(3.39) 



and it has already appeared in (|2.20|) . The gluon coefficient functions can be expressed 
in terms of <M as 



* 9 n(6) = K-x(6, 6) + $Li(6, -6 - 6)] + (6 - -£ 
2 \ 2(n- 1 



(3.40) 



<^(6) + 



n(n-2) y nv "' y ' n(n-2) 
Explicit expressions for a few first moments read: 



K_i(-6-6,6) + ^- 1 (6 + 6,-6 



$1 



fi§ ff = 



$1 = 5(6-6 



Of 



31 



(6-6), 



(3.41) 



*? = 14 (6 3 - 6 3 ) , = ~ (59 6 3 + 66 2 6 - 666 2 - 59 6 3 ) ■ 



(3.42) 



The expressions (|3.26|) , (|3.30| ) and Q3.37| ) for the one-loop gluon contribution to the 
antisymmetric part of the T-product of two electromagnetic currents ( |2.3|) and the struc- 
ture function g2(xB,Q 2 ) present the main result of this section. 



4. Properties of the twist-3 contribution 



According to ( |3.37| ), the moments of the structure function are given by integrals of 
the quark-gluon and three-gluon distribution functions, -D g (6) and D g (£i), respectively, 
over momentum fractions of partons with the weights defined by the coefficient functions 

^te),^te)and^fe). 

The quark coefficient function (|3.39|) vanishes at n = 1 and for higher moments 
is a homogenous polynomial in momentum fractions of degree n — 3. For the gluon 
coefficient functions, Eqs. Q3.4Q) , one finds that f^=i(£i) = 3>n=i(&) = and, therefore, 
the gluon contribution to the first moment of g2{x,Q 2 ) vanishes, in agreement with 
the BC sum rule, provided that the three-gluon distribution function D g (£i) does not 
have additional singularities. Moreover, in contrast with the quark coefficient function, 
f^_ 3 (£j) = $^ =3 (£j) = 0, and gluon contribute to the structure function starting from 
n = 5th moment. For n > 5 the coefficient functions and $^ are homogenous 
polynomials in £j of degree n — 4. 

The explicit expressions for the lowest moments Q3.37| ) look as follows 



q 



(4.1) 
(4.2) 



x 



(3^ - 2^6 + &)D q {& Q 2 ) - - 6)^(6; Q 2 ' 



7T 1440 



(4.3) 



x 



[(5£ - 4£ 3 £ 3 + 3£ 2 £ 3 2 - 26^3 + €t)D q {& Q 2 



7T 



3091. - 207 ,~ _ 2 , 



2240 



1120 



D 9 (^Q 2 



Using (|2.18f) and ( |3.36|) we express the moments in terms of reduced quark and gluon 
matrix elements 



dxx 2 gf-\x,Q 2 ) = -^eJ((5 °(Q 2 ))), 



dxx A g^~ 3 {x,Q 2 
2 



(4.4) 
(4.5) 



1 

[ dxx 6 gl w - 3 (x,Q 
Jo 



rv 847 

((S° 2 (Q 2 ))) - 2((S l 2 (Q 2 ))) + 3((S 2 (Q 2 ))) + ^^((G° 2 (Q 2 ))) 



(4.6) 



7 E e S [(( 5 4 (Q 2 ))) - 2((5](g 2 ))) + 3<(s 4 2 (q 2 )» - 4((s 4 3 (g 2 ))) + 5((5 4 4 (g 2 ))) 



-1 .-7 



+-( 3 ^((Gl(Q 2 ))) + 2 ^((G\(Q 2 ))) 



71 



560 



We would like to stress that the relations (13.371 ), (|4.1|) and (|4.4j ) take into account 
the leading order contribution of the three-gluon operators and they do not include 
0(a s )— corrections to the coefficient functions of quark-antiquark-gluon operators. The 
latter corrections have been recently calculated in |33| using a different operator basis. 

We can make one further step and define the momentum fraction distribution of the 
transverse spin carried by gluons in the nucleon by the expression similar to Q2.11| ): 



dv(p,s\ [u + v)0^{y , u, — u) + 2u 0^{u, v , — u) + {u — v)Op(u, — u, v ) \p, s) ■ 

i 

= 2 [s,{px)-p,{sx)]j\e 2 ^ (■ Ag T {i^ 2 ) 



(4.7) 



so that (cf. {WM ) 



i j -i 



(4.8) 



and Agx(^) = Ag^(— £). It is easy to see that the first contribution in the square brackets 
in (|3.37| ) can be easily rewritten in terms of AgT(n,Q 2 ) = J 1 (i^ n_1 A(7y(^). Unfortu- 
nately, the two coefficient functions $^(£i) and f^ 9 (£j) do not coincide and, therefore, the 
full gluon contribution to Q3.37|) cannot be, strictly speaking, reduced to the contribution 
of Agx{C,) alone, as a yet another manifestation of the fact that we are dealing with a 
three-particle problem. 

Our main observation is that such a reduction can, nevertheless, provide a reasonable 
approximation to the moments g2(n, Q 2 ) 



L 



W A,(&; Q 2 ) CW*) « c(n) / T>{D g (& Q 2 ) , 
i J-i 



(4.9) 



with the coefficient of proportionality c(n) that is independent on the momentum transfer 
Q 2 . In this case, moments of the structure function (|3.37p can be expressed in terms of 
the quark and gluon distributions as 



j dx B x n B 1 g 2 {x B , Q 2 ) w ^ e l J dx 



(4.10) 



x [Aq T (x B , Q 2 ) + Aq T (-x B , Q 2 ) + Cf ~\n) Ag T (x B , Q' 

where Aq T and Ag T were defined in (|2.11| ) and (|4.8|) , respectively, and 

a s 1 



C* w " 3 (n) -- 
9 y 1 2nn+l 



[(V>(n) + 1e + 1) c(n) - 1]. 



(4.11) 



The reason for ( )4.9|) to hold is that, as we shall argue below, the both projections of the 
three-gluon distribution defined in the l.h.s and the r.h.s of Eq. (4.9) can be identified 
to a good numerical accuracy with the contribution of the three-gluon multiplicatively 



renormalizable operator with the lowest anomalous dimension. To be more precise, this 
statement refers to the "purely gluonic" operator, defined without taking into account the 
mixing with the quark- ant iquark-gluon sector. We have checked that the mixing between 
three-gluon and quark-antiquark-gluon operators does not modify the result significantly; 
a detailed analysis will be presented elsewhere [j34 |. 



To justify this, we note that in a "purely gluonic" sector an arbitrary multiplicatively 
renormalizable twist-3 three-gluon operator O^^ can be characterized by the coefficients 
in its expansion over the basis of operators [G>]jv defined in ( |3.34|) 



[N/2]-l [JV/2]-l 

N , a = £ ™ k N, a [G»] k N = - £ < Q (tt-PX" 1 "*) ( 4 - 12 ) 



k=0 k=0 

or, equivalently, by a characteristic polynomial 

[N/2]-l 



6) = \ £ <« (tf^- 1 -* - ef- 1 -^) • (4.i3) 



k=0 



The subscript < a < [N/2] — l enumerates the operators and we assume, for definiteness, 
that the operators On, a are ordered in such a way that a smaller a corresponds to a lower 
anomalous dimension. With this definition, it follows from (|3.36 ) that the reduced matrix 



element of a multiplicatively renormalizable operator On,o is given by a weighted integral 
of the three-particle gluon distribution: 

((0 N , a )) = j W W N , a {£i) D g &) . (4.14) 

For the purpose of our discussion, we assume that the expansion coefficients w Na in 
( |4.12| ) are calculated by an explicit diagonalization of the mixing matrix given in || so 
that the polynomials WN,a{£i) are known functions. 

To prove our assertion, we have to show that the coefficient functions and 
Q^ 9 ((,i) are numerically close to VFjv )a= o(£i), a ^ least for sufficiently large values of n — 
N + 3, or, equivalently, the both sides of (|4.9|) receive a dominant contribution from 



({On,oc=o))- This is not straightforward since the characteristic polynomials Wjv,a(£i) for 
different a are not mutually orthogonal with respect to any simple weight function, the 
reason being that the mixing matrices are not symmetric. 

In order to make a meaningful comparison we use the conformal symmetry that allows 
to rewrite the mixing matrices in a different basis such that they become hermitian (see 
55 1 for details). In the present context, the idea is that the conformal symmetry 



allows for a unique analytic continuation of the functions Wjv,a(£i)) ^jvte) anc ^ 
defined on the plane £i +£2 + £3 = to the full three-dimensional space of the momentum 
fraction variables £j 

WV,«(6,60 => Wiv, Q (6,6,&) ( 4 - 15 ) 

etc., such that Wjv,a(£i, £3) = Wjy,a(£i) f° r £1 + ^2 + ^3 = and the functions Wjv,a(6i) cor- 
responding to different multiplicatively renormalizable operators are mutually orthogonal 
with respect to the conformal scalar product 

(W n JW Ni(3 ) = I [d£\ W N ^ t )W N ^) ~ 8 aP , (4.16) 
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Figure 3: Conformal projection (see text) of the coefficient functions <3?n(£i) (crosses) and 
^n 9 (£i) (diamonds) on the contributions of multiplicately renormalizable three-gluon operators 
enumerated by an integer variable a in the order of the increasing anomalous dimension. Mixing 
with quark-antiquark-gluon operators is neglected. The two plots correspond to the moments 
n = 13 (right) and n = 33 (left), respectively. 



where the integration measure is defined as 

[d£] = «2^ 3 5(6 + 6 + 6 - 1) • (4.17) 

Using this method, we can expand and f^- 9 (£j) over the set of orthogonal polyno- 

mials Wjv,o:(£i) to arrive after reduction to £i + £2 + £3 = at 

[JV/2]_1 i$s\w \ [JV/2M m 

*«&) = E = ll*ill E ^Srr ^ 4 - 18 ) 

a=0 II WjV,a|| a =0 II^A^H 

with n = iV + 3 and similarly for f2* 9 (£j). Here, the norm is denned in a natural way as 
The results of our calculations of the (normalized) expansion coefficients 

(h\w N , a ) wfe> _ (419) 



W&4 • \\W N , a \\ • \\W N , a \\ 

corresponding to the functions Q 9 n and Q^, respectively, are shown in Fig. |3| for n — 13 
and n = 33. It is clearly seen that the sum ( |4.18|) is dominated by the single contribution 
a = of the three-gluon operator with the lowest anomalous dimension. The quality of 
such approximation is improving for larger moments n. Moreover, we find that, in accord 
with (|4.9|), the two functions and f^ s (£j) are, in fact, close to each other with the 



coefficients c(n) in ( |4.yQ and ( J4.11|) given by 



(*S|H'k.o) 



Since the contribution of f2^ s (£j) to ( |3.37|) reflects the mixing of three-gluon with quark- 
antiquark gluon operators, this is an indication that the observed dominance of the lowest 
three-gluon state is not obstructed by this mixing. 
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Figure 4: Conformal projection c(n) of the full gluon contribution, Eq. ( 4.2C| ), onto the trans- 
verse spin gluon distribution defined in (|4.7| ). 



From Eq. ( |3.40|) it is easy to see that c(n) — > 1 for large moments n whereas from 
( |3.41| ) it follows that c(5) = 31/25. The values of c(n) for other (odd) moments n < 25 are 
shown in Fig. The n-dependence is very smooth (for odd n) and can be approximated 

as 



c(n) = 1 + 11.11— -25.74 — . 



(4.21) 



Analytic expressions for c(n) can be worked out in the large- N c limit and will be presented 
in M. 



5. Conclusions 

We have presented a detailed calculation of the three-gluon twist-3 contribution to the 
flavor-singlet structure function gz(x, Q 2 ) to the one-loop accuracy. The result is encour- 
aging as it indicates that the gluon coefficient function is close to that of the three-gluon 
operator with the lowest anomalous dimension, at least for large moments. This allows to 
hope for a simplified description of the scale dependence of g?,{x, Q 2 ) in terms of DGLAP 
equations, similar as for the structure functions of leading twist. The construction of 
this approximation requires a more detailed analysis of the evolution equations for the 
flavor-singlet twist-3 operators and will be given in a forthcoming publication f34| . 
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Note added 



When this work was in preparation, the work [36| appeared with the calculation of 
the singlet twist-3 coefficient functions using a different approach. The relation of this 



calculation to our result is not obvious because of different operator basis. It appears 
that the answer for the n = 5 moment of #2 0*0 given in (|4.1|) is in agreement with the 



appropriate projection of the coefficient function calculated in [ 36fl . We thank A. Belitsky 
for the correspondence on this topic. 



A Appendix 



The short-distance expansion of the nonlocal two-gluon operator looks like [] 

oo ^ 
n=l ^ '' 

The twist- 2 contribution to the r.h.s. originates from the local operators completely sym- 
metric and traceless with respect to their Lorentz indices 







t n _ \\u n _|_ i) q x i x vi x t l 2 ■ ■ ■ x nn x vn+i) Gfj, ia D^ 2 ■ ■ • D/i„G Mn+iaf (0) , (A. 2) 

n=l ' v 



where Dp = Dp — Dp- It is straightforward to verity that the same expression can be 
obtained through the following integral representation 



[G xa (x)G va {-x)\ 



1 d 

!v 1 ' duu— — G xa (ux)G xa (— ux) . 
o ox i, 



(A.3) 



Subtracting the twist- 2 contribution from the nonlocal two-gluon operator we obtain the 
relation 

G xa (x)G va (-x) - / udu- — G xa (ux)G xa (-ux) = 



udux r 



dx v 




d 

G xa (ux)G va (-ux) - - — G xa (ux)G pa (-ux) 



OXp ox v 



+ / udu 
'o 



G xa {ux)G ua (-ux) - G ua (ux)G xa (-ux) 



(A.4) 



which can be expressed in terms of three-gluon operators using the QCD equations of 
motion as follows. 

We begin with the expression in the second line of ( A.4 ): 



d d 

Gp a (ux)G ua (-ux) - —G ua (ux)G pa (-ux) 



dx 



G 



pa 



uDp — uDp — i I vdvG xp (vx) 



G va (-ux) - (p u) . (A. 5) 



7 Throughout this Appendix we assume that gauge factors in the adjoint representation are inserted 
in between the gluon fields. 



on 



Following || , we introduce a derivative over the total translation 

d r 

d p G^ a (ux)G ua (-ux) = — G^iux + y)G va (-ux + y) 



dy, 



G^UX) 



Dp + Dp-i I dvG X0 (vx 



' xp \ 



G ua {-ux) (A.6) 



so that 







G m (ux)G ua (-ux) - ud p Gp a (ux)G ua (-ux 



dx p ~» a 



(p <-> u) 



= —2uG xa (ux)D a G ux (—ux)—i / (u — v) dv G xa (ux)G xl/ (vx)G xa (— ux) 

J —u 

= 2ud a G xa {ux)G xp { — ux) — 2uG xa {ux)D a G xv {— ux) + 2iu / dv G xa (ux)G xa (vx)G xu (—ux) 

J —u 

/u 
(u - v)dv G xa (ux)G xu (vx)G xa (-ux) (A. 7) 

■u 

where we used the Bianchi identity DjG pa — DpG ua = DcPup to arrive at the expression 
in the second line. Neglecting total derivatives and terms ~ D a G au this gives 



udu x, 



d d 

G xa (ux)G ua (-ux) - - — G xa (ux)G pa (-ux) 



dx 



dx y 



(A.f 



i / udu / dv 



2uG xa (ux)G xa (vx)G xu (-ux) + (v - u)G xa (ux)G xu (vx)G xa (-ux) 



Next, we consider the antisymmetric combination of gluon fields in the last line in 
( |A.4j ) and rewrite it as 



G lwt (x)G va (-x) - (fJL<-+ v) 



1 du d 



U OX\ 



Gp a (ux)G va (-ux) -(/jh-+v) 



duXy 



Gp a (ux)D\G ua (-ux) - Gp a (ux)D\G ua (-ux) -(//<->■ i/) . (A. 9) 



By a repeated application of the the Bianki identity and separating (and then neglecting) 
contributions of total translation and equation of motion terms ~ D a G au , one arrives 
after some algebra at the following expression: 

-1 r u 



x. 



G m {x)G va {-x) -(n++v) 



2i du dv 

'0 J-u 



G xa (ux)G xa (vx)G xl/ (-ux) 



- G xa (ux)G xu (vx)G xa (-ux) + G xl/ (ux)G xa (vx)G xa (-ux){A.10) 



that yields 



udu 



G xa {ux)G va (-ux) - G ua (ux)G xa (-ux) 



1 

2 



du (1 — u ) I dv G xa (ux)G xa (vx)G xu (—ux) — G xa (ux)G xu (vx)G xa (— ux) 

+ G xu (ux)G xa (vx)G xa (-ux) . (A.ll) 



Taking the sum of the expressions in ( A.8 ), (|A.11|) and the color trace, we obtain the 
result (|3.11|) quoted in the text. 
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